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INTRODUCCIÓN 


Normalmente hablando, el aire natural, además de estar formado por 
una mezcla de diversos gases en diferentes concentraciones bien 
conocidas, contiene una parte constituida: por partículas en suspen- 
sión (aerosoles) de diferentes tipos, naturaleza y proporciones. Por 
ello, el aire natural perfectamente limpio, en términos reales, no 
existe y por el contrario, como algo común, siempre está sucio o 
tiene algo de suciedad que puede ser de origen natural o antro- 
pogénico; dependiendo de las circunstancias, este sucio puede ser 
inofensivo, benigno y hasta positivo, o puede ser dañino y maligno. 
En el sentido positivo, existen procesos convenientes en la 
naturaleza (ciclos) que ocurren en diferentes partes de la geoesfera, 
que dependen en gran parte de la circulación e intercambio de 
partículas naturales y gases en la atmóstera. Así la calificación y 
cuantificación de la suciedad del aire por partículas y/o gases 
extraños o no, queda sujeta a su origen, papel, nivel, alcance, 
efectos y consecuencias que de ella se deriven. Se puede hablar, 
-entonces, de la suciedad combinada natural (inofensiva o dañina) y 
de la suciedad artificial (generalmente dañina) del aire. Cuando está 
suciedad, combinada o no, se presenta en el aire por encima de 
ciertos niveles o cantidades, o está compuesta por elementos 
extraños (gases o partículas) a su composición normal, podremos 
hablar de una contaminación general del aire en el sentido negativo 
de la expresión. o 


Específicamente, el aire está contaminado cuando la 
concentración de cualquiera de sus “ingredientes secundarios 
(naturales y/o artificiales) es lo suficientemente elevada como para: : 
(1) dañar la salud de las personas y otros seres vivos, (2) causar 
pérdida de bienestar, (3) disminuir la visibilidad y/o (4) producir 
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cambios climáticos (OPS, 1980). La pérdida de bienestar puede ser 
por: (a) acumulación de polvo y suciedad, (b) emanación de olores 
desagradables y/o (c) disminución de los rayos solares. 


La contaminación del aire como problema ambiental central 
abarca una serie de posibles escenarios, formas o manifestaciones 
que varían tanto en dimensión espacial como en dimensión temporal. 
Se reconoce, entonces, que este problema no es uno solo sino 
muchos interconectados, con diferentes características que lo hace, 
al momento de su estudio, conocimiento y solución, un problema muy 
complejo de naturaleza inter y/o multidisciplinaria. 


Una visión general de las diferentes categorías y escalas 
espaciales y temporales de cómo este problema se presenta, 
muestra que la contaminación del aire puede ser un episodio o 
evento muy localizado y rápido, pero muy agudo o grave, o llegar a 
ser un fenómeno progresivo que potencialmente puede estar 
afectando globalmente a todo el planeta. 


Al respecto, la Tabla | muestra como la contaminación del aire 
puede ser categorizada tomando en cuenta cuatro aspectos gene- 
rales: la escala horizontal, la escala vertical, la escala temporal y la 
escala de pertinencia o competencia jurisdiccional. 


"HORIZONTAL" "VERTICAL - “TEMPORAL - JURISDICCIÓN 
Local ., < Altura de Chimeneas Horas Municipal-Vecinal 
Ubana * -1km Días Distrital 
Regional _Troposfera Meses . Estatal-Nacional 
Continental Estratosfera Años Internacional-Nacional 
Global: -Atmósfera Décadas Intemacional 


TABLA 1. Ámbitos bajo los cuales la contaminación del aire puede ser considerada 
(adaptado de Stern et al., 1984). 
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La primera escala (horizontal) indica cuánta superficie de la 
Tierra está envuelta en el problema. La segunda (vertical) indica qué 
profundidad (o altura) de la atmósfera está involucrada en' el 
problema. La tercera (temporal) señala cuál puede ser la duración del 
problema y eltiempo disponible para que se resuelva. La cuarta y 
última, determina cuál sería el nivel de los organismos e instituciones 
jurisdiccionales a los sión le compete atacar y resolver el problema 
oficialmente. 


La focalización del problema de acuerdo a este esquema 
define, entre una gran variedad de cosas, su evaluación teórica y/o 
práctica, pronósticos o predicciones, recursos, costos, impacto o 
consecuencias, estrategias (de mitigación o adaptación), legislación, 
etc., que evidencia una gran complejidad. 


Hoy día los ejemplos de contaminación . ambientál (ds < origen 
antropogénico), a diferentes escalas y dimensiones, y en particular 
los que tienen que ver con el aire y que lo someten a traumas 
adicionales, además de los naturales, abundan por doquier debido al 
alto grado de desarrollo tecnológico que ha alcanzado la civilización 
modema contemporánea. Por mencionar algunos, que destacan por 
su relevancia y trascendencia, se pueden. citar el accidente de la 
planta nuclear de Chernobyl, Ucrania, (contaminación radioactiva 
letal), lluvia ácida, episodios de altísimo grado de contaminación 
(smog fotoquímico) en grandes ciudades sobre-pobladas, accidente 
industrial de Bophal, India, (contaminación química venenosa), 
destrucción de la capa estratosférica de ozono o de la ozonósfera 
(penetración de radiación solar ultravioleta en exceso), el 
recalentamiento global (aumento de la temperatura media global). 


De los anteriores, el problema conocido ampliamente como 
“Efecto Invemadero” o “Recalentamiento Global” (en inglés 
“Greenhouse Effect”) y el problema de la contaminación 
antropogénica del aire de una gran ciudad, serán en este curso 
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«discutidos a un nivel básico para introducir conocimientos. lid 
tales pen lleven al entendimiento esencial de los mismos.: 


El primero de estos problemas, contarme « a la Tabla L abarca 
a: 106 el planeta y a toda la: atmósfera, es a. largo plazo - -y de 
pertinencia internacional, es decir, nos concierne a todos. La historia 
:de.cómo este problema llegó a ser identificado y cuantificado, en sus 
primeras. aproximaciones, puede ser encontrada en. Jones 8 
Henderson-Sellers (1990) y RSAS (1997). Una discusión sobre el 
tratamiento y. «cobertura . dado al fenómeno” por. los. medios de 
-comtinicación “social: impresos anglosajones entre 1987 y 1990,. 'es 
«presentada: por Wilkins-:(1 993). Y el tema y sus consecuencias es - 
tratado a' un buen nivel general por Aubrecht (1988), Sagan (1990), 
Mitcheil:(1990), Erickson (1992) y Alonso 8 Ramis (1996). Un trabajo 
más completo que abarca aspectos adicionales sobre 2 cambio 
: dlarcicn do es el de Vázquez-Abeledo q 998). 





Eb sudo problema, de “acuerdo a a misma tablas Abarca 
A Er una zona pequeña de tipo urbana hasta 1 km de altura en la 
atmósfera, puede durar algunos días y le. compete jurisdiccional- 
mente: a. la“instancia distrital. Seguramente este problema hizo su 
aparición en'las primeras décadas. «de la era industrial durante el siglo 
XIX; cuando a' partir de 1860 la calidad del aire urbano de algunas 
ciudades : europeas “comenzó a ser “sometido a. drásticos cambios 
nunca: antes vistos. Por 1911, Aitken (citado por-Davies, 1987) ya * 
había: mostrado Ja preocupación sobre las enormes cantidades de 
sulftro:: inyectadas' a la atmósfera, debido..a.la quema de carbón 
natural: para uso doméstico e industrial. Hoy día. la contaminación 
atmosférica: de las* -grandes ciudades del.. mundo. .es un problema 6 
común. «Aun cuando ambos problemas comenzaron a ser yistos 


seriamente a partir del siglo XIX, el segundo es.mucho más viejo que 
el api (Halliday, 1961). 
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: Para entender el recalentamiento global terrestre por causas 
antropogénicas, el cual se presume está “en progreso, 
necesariamente hay que entender primero el efecto de invernadero * 
natural terrestre que hace a este planeta habitable y confortable para 
la vida, en .comparación.con los otros planetas terrestres,. en 
- particular, con Marte y Venus. Esto se hará con una matemática y 
una física muy simples. no Ea 
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EL EFECTO DE INVERNADERO NATURAL TERRESTRE 


¿Es obvio que la Tierra como cuerpo en el espacio, en promedio, tiene 
_Qque estar en equilibrio de radiación (equilibrio radiativo); de lo 
contrario, su temperatura iría en aumento hasta que explotara o se 
evaporara, o su temperatura iría disminuyendo hasta que se con- 
gelara. 


En tal equilibrio se cumple que: 
AF=FÍ- FI =0, (1) 


donde Fl es el flujo de energía solar que recibe la tierra y que entra 
por su atmósfera y FT es el flujo de energía térmica terrestre que 
sale irradiado hacia el espacio exterior. 


La Tierra, a una distancia promedio del Sol, recibe de él la 
cantidad de energía S, = 1368 joules/(s m”) = 1368 Wm”, la cual es 
intersectada por la cara de la Tierra que está dando en ese momento 
al Sol; esa cara (disco) tiene un área de A, donde R, es el radio de 
la Tierra. Ahora bien, la Tierra refleja parte (fracción del total) de esa 
energía y la devuelve al espacio. Esto se conoce como albedo, a, el 
cual depende de las propiedades del sistema atmósftera-superficie. 
Entonces la parte de S, que entra a la tierra es (1-a), y el total que 
entra a través de toda la cara es, 


FS =8S(1-— aJrR?. (2) 


Al mismo tiempo, la Tierra emite al espacio exterior (como un 
cuerpo estérico calentado a una temperatura radiación, T.), una 
radiación térmica, por segundo y por metro cuadrado, dada por la'ley 
de Stefan-Boltzmann (ver más abajo): oT/, donde o es una 
constante igual a 5,6705 x 10” W/(m' k). esta radiación es emitida a 
través de la superficie de todo el globo terráqueo que es igual a 
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Ar R;. por lo tanto, la cantidad se AMEnalA traLclkacla al espacio por 
toda la Tierra es, co 


y FI = 4x1 R207.. (3) 
Por condición de la ec. (1): 
8.(1 - añ; =41Af 07, (4). 


la cual se reduce a, 


S.(1- aa =oT,. (5) 


Despejando 7, de la ec.(5), 
T, = [(1/40)8,(1 - aJ]"".(6). 
con a= 0,35 y demás valores dados arriba, T, = 255 K (18 20). 


- Este valor de T, representa la temperatura promedio corres- 
pondiente a la energía térmica emitida tanto por la superficie de la 
Tierra como por su atmóstera, vista desde el espacio. Esta radiación 
se detecta como si viniera desde una altura alrededor de los 5 km 
sobre el nivel del mar. Así el nivel de radiación efectiva del sistema 
combinado no se localiza en la superficie sino ácida en 
la mitad pp Sea de la atmósfera.  : 








Con la ec. (5) podemos cuantificar el: lujo da energía global 
anual involucrado en el equilibrio o balance radiativo entre la Tierra y 
el espacio. Así, con el mismo valor de a y S,,. FT = 222,3 Wm”. La 
energía solar que llega a la Tierra, a una distancia promedio del Sol, 
es, como se ha indicado antes, S, = 1368 Wm?. Pero de la ec. (5) se 
ve que sólo S,/4 (= 342 Wm”) es el valor que entra en el cálculo del 
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- balance, debido al factor (1/4) el cual proviene del cociente entre la 
sección transversal de la esfera terrestre (RA?) y su superficie 
(4xR?). Este factor geométrico es fácilmente interpretado si 
consideramos la rotación de la Tierra. En promedio anual, un lugar 
sobre nuestro planeta es iluminado por 12 horas y no por 24; 
además, está iluminación diaria no es siempre perpendicular, sino 
incide de acuerdo a la posición del Sol en el cielo. En promedio, el 
coseno de este ángulo de incidencia es de aproximadamente Y, y 
como la iluminación es de medio día, o sea, %, el flujo recibido por la 
Tierra es en promedio global anual de 342 Wm” (S,/4). En el caso de 
a = 0 (se absorbe todo y se emite todo), en balance radiativo, Fl = 
Fl = 342 Wim” (Fig.1), pero como a = 0,35, FT =F] =-222 Wm”. La 
diferencia de 120 Wm* (aS/4) corresponde al flujo que es devuelto al 
espacio exterior por reflexión (Fig.2). En estos dos ejemplos hemos 
considerado que la atmósfera terrestre no interviene (o no existe) en 
el proceso, por lo que el albedo corresponde a la superficie. 





FIGURA1 |  FIGURA2 


Supongamos 'ahora un. tercer escenario donde se absorbe 
todo (a =.0), pero ya la atmósfera interviene en el proceso de balance 
tadiativo. Bajo esta condición, la atmósfera emite 342 Wim” hacia el 
espacio exterior, pero también hacia abajo (Fig.3). Y con los. 342 
WmY que recibe la superficie, ésta última emite hacia arriba 684 Wm?* 
(Fig.3). En un cuarto escenario suponemos lo mismo, excepto que 





tomamos en cuenta los 120 Wm” reflejados al espacio; en éste caso, 
lo qu > llega a la superficie es 222 Wm” y lo que debe salir al exterior 
(desde la atmósfera). es lo mismo por condición de: equilibrio radiativo 
(Fig 4). Por otra parte, la atmósfera devuelve, en equilibrio, también 
222 Wm”, que. son recibidos por la superficie junto con los 222 Wm” 
qué “entran desde el exterior, para ún total de 444. wm? que. son 
emitidos desde la superficie 1, 4). "e 











- FIGURA3. Aquí.T, =7, y T,= 7, 
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La pregunta que sigue es: ¿cuál es entonces la temperatura 
de radiación, T., de la superficie terrestre y cuál es la Isern de 
radiación, T,, de la atmósfera? 


La atmósfera, por encima del 99%, está compuesta de 
moléculas diatómicas y triatómicas las cuales, absorben radiación 
(que afecta sus modos de vibración y rotación) en un amplio rango 
del espectro electromagnético. En particular, ellas son más efectivas 
absorbiendo en el rango infrarrojo, IR, (A = 1-100 um) que en el 
visible y en el cercano IR (A = 0,3-1 jm). Como resultado, el grueso 
de la radiación solar (la cual está concentrada alrededor de este 
último intervalo), pasa a través de la atmósfera y es absorbida por la 
superficie, la cual'se calienta hasta alcanzar, en promedio, una 
temperatura de equilibrio 7, 


En el equilibrio, la superficie irradia hacia arriba, hacia la 
atmósfera, energía IR la cual es parcialmente absorbida por estas 
moléculas de la atmósfera como el H,O, CO,, O,, etc. (gases 
invernadero). En equilibrio radiativo, esta energía parcialmente 
absorbida es re-emitida hacia arriba y hacia abajo (hacia la 
superficie). Esta absorción es selectiva en lo que se refiere a A 
(longitud de onda de la radiación), pero por el momento supongamos 
que la misma no lo es, por lo que ahora se considera que la 
atmósfera absorbe como un cuerpo “gris” que tiene una emisividad e, 
<1,e independiente de 4 (en el IR). La emisividad es la fracción de 
energía emitida (por unidad de tiempo y por unidad de superficie) por 
el cuerpo con relación a la que emitiría si irradiara en todas las 
longitudes de onda (cuerpo negro). Sie, = 1, entonces el cuerpo 
emite como “cuerpo negro”. (en todas las longitudes de onda). La 
distribución de esta energía en función de 4 genera una curva la cual 
es descrita matemáticamente por una función que al integrarse (área 
bajo la curva), da la "energía total emitida en todas las longitudes de 
onda (por unidad de tiempo y de superficie). Esta última, es- 
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proporcional a la cuarta potencia de la temperatura de radiación, 
llamada temperatura efectiva 7, (ley de Stefan-Boltzmanmn: o T;' ). 
Con estas aproximaciones y consideraciones, se puede 


construir un modelo del sistema atmósfera-superficie de una sola 
capa, como se muestra en la siguiente Figura 5 (Ackerman, 1991). 


| (S/4N1 a) | 


A 

A ATMÓSFERA 

py | 

y Bl oT! 
SUPERFICIE 

FIGURA 5 


La atmósfera es representada por una sola capa con una 
temperatura uniforme 7,. En el equilibrio, el balance de energía en 
este modelo simple puede ser representado por un sistema sencillo 
de dos ecuaciones con dos incógnitas: 7, y 7.. 

En la superficie: 

(S/4)(1 - a) + e,0T, =0T;. (7) 


En la atmósfera: 


£,oT/ + e,0T/! =€,0T.*. (8) 
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Resolviendo el ió 
o T, =[11(2)"17,, (9) | 
| T, = [(1/498,1 — aylo(t — 1/26,)1=T,/[1— (1/2), 1%, (10) 
T,=TA2(1 - 1/26,)1 = T/1(2— EJ”. (11 , e 
Las ecs, (10) y (11) dan, SEpccitvamente, T, y T, en función de : 


T, y de e, con 0 < e, < 1. La Fig. 6 muestra gráficamente la variación 
- de 7, y T, en función de e, (con T, = 255 K). 
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FIGURA 6. Temperaturas atmosférica y superficial en función de la emisividad 
atmosférica (modelo de una sola capa) 
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Note-que Y, 5.4 que T.S Fo 


* La flecha inferior. en la Fig.6 Aka 1991): indica: el valor 
des 7 “que es la emisividad de la atmósfera espectralmente 
_promediada. Los correspondientes 'valores * de “Ty T, son, 
respectivamente, 288 K (flecha SUpSñtN y: 244 K (lscha intermedia), 


“También de la Fig. 6 se obsetvá que para:s, 0, T, qa Ts 
255 k y T, >214 K; parae¿>1, T.> 303 K y .T, > T, = 255 K. 


Para £, = 0 (no hay atmósfera), el valor: exacto de T, no tiene 
significado físico. : o: : ty 


Actualmente la temperatura de la. superficie terrestre está 33 
K (= 288 K -— 255 K) por encima de la temperatura de equilibrio 7, y la 
temperatura de la atmóstera está 11 K por debajo (= 255 K — 244 K). 
El calentamiento de la superficie se debe al flujo que devuelve la 
atmósfera hacia abajo. Esto último sé conoce con: el nombre de 
calentamiento global o efecto de invernadero natural terrestre. j- 


Sin este efecto, la superficie de nuestro planeta sería signifi- 


cativamente más fría y el ambiente menos hospitalario y menos 
apropiado para la vida. . 


Problema 1: 


Calcular T, y T, para los casos planteados en las Figs. 1-4. 
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ALTERACIÓN DEL EFECTO DE INVERNADERO 
(CALENTAMIENTO GLOBAL ADICIONAL) 


Las concentraciones de CO, y H,O determinan fundamentalmente, 
junto otros gases-traza que absorben en el IR, la emisividad de la 
atmósfera terrestre. Si todas estas concentraciones o cualesquiera 
de ellas se altera, la temperatura de la superficie se alterara como 
también lo hará la temperatura atmosférica (ecs. 10 y 11). j 


Para ver cómo se altera la temperatura de la superficie, 
rearreglemos la ec. (10) como, 


[1 41/2),Jo7 = (1/4)5,(1 - a). (12) 


- Diferenciemos esta expresión con respecto a la temperatura super- 
ficial y a la emisividad: 


[1 1/2)€, 40 T AT, + 07,(-1/2)€,Ae, = 0. (13) 

Busquemos cuánto varía la temperatura superficial, AT, en 
función de un cambio As,. dividiendo por T,; ambos miembros y 
rearreglando, 

[1 —(1/2)e, JAS(AT/T,) - 0(1/2)e,As, =0 > 
AT, = T, £,Ag, / (8 — 4£,). (14) 


La Figura 7 muestra gráficamente la ec.(14) para diferente 
valores de e,. 
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FIGURA 7. Cambio en la temperatura superficial terrestre en función del cambio de 
- emisividad. 


Se estima (Ramanathan y Coakley, 1979) que si se dobla en 
el futuro la concentración actual de CO, (-358 ppmv > 2x 358 
ppmv), 


producto de la actividad antropogénica, habrá un flujo adicional 
“atmosférico hacia debajo y a nivel de la superficie de 4 W/m', 
equivalente a un aumento de la emisividad de 0,025. La relación 
entre la emisividad del CO, y su concentración en la atmósfera [CO,] 
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ha sido dada por Staley y Jurica (1972), con referencia al nivel de los 
1000 mb de presión, quienes han establecido la relación e, = 
0 ,0235LnN[CO, ] + 0,0537: Aplicando esta relación a [CO,] = 358 ppmv 
y [CO,] = 716 ppmv, Ae, = 0,020. Del gráfico anterior, se ve que esto 
equivale a un cambio de la temperatura superficial terrestre de AT, = 
1,2 K.. Una pregunta que inmediatamente surge es, ¿Para cuándo se 
alcanzará el doble de la concentración de CO, en la atmósfera? Esto 
es materia de extenso e intenso estudio a través de modelos : 
matemáticos. 


Problema 2: 


Encontrár cómo varía AT, VS Ae, y comparar con la curva dada en 
la Fig.7. 


. La predicción hecha por este modelo sobre-simplificado 
se ba basado, entre otras, en las siguientes aproximaciones: 


hi La atmósfera tiene una sola capa homogénea. 
2. - Toda la energía solar llega a la superficie. 
3. La atmósfera emite como un cuerpo gris. 
pa La capa contiene sólo CO,. 
-5.-La atmóstera no contiene nubes. 
6.- La radiación es el único mecanismo de transferencia de energía. 


7.- La radiación solar cae perpendicularmente. 
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8.- La superficie absorbe toda la energía que le llega. 


En general, una variación de la temperatura T, y de T., 
significa una variación también de 7,, la cual llamaremos AT,. Para 
estudiar, bajo un enfoque alternativo, esta variación debido a una 
alteración de la composición de la atmósfera por contaminación, 
consideremos que el globo terráqueo tiene una temperatura de 
radiación 7,, asociada con un flujo hacia el espacio exterior dado por 
la ley de Stefan-Boltzmamn (S, = 0T,). cuando el globo terráqueo 
adquiere una nueva temperatura de equilibrio 7,, el flujo ha variado 
de S, a S,, en una cantidad AS, = S, - S,. En símbolos, 

S,>8,>AS,=8, -S,>8,=AS, +8, (15) 
Como consecuencia, 
T,> T,>AT,=T, - T,> T, =AT, + T,. (16) 
De ahí que: 
AS, > AT. (17) 


¿Cómo se relacionan estos dos cambios? Siguiendo a Knox 
(1999), aplicamos la Ley Stefan-Boltzmann a S,, y por aplicación de 
ec.(15) y (16). 

S, =07,* > (S, + AS,) = 0(T, +aTy, M8)" 
Dividiendo ec.(18) por S,, 
[1 +(4S/S)] =(0/S XT, + AT)", 
[1 +6] =(1/7,) (T, + AT”, (19) 


Donde e = (AS//S,) representa la variación de S, como fracción del 


22 | Introducción a la Física de la Contaminación Atmosférica y del Cambio Climático 


flujo original S,- Despejando AT, en ec.(19), 
[res (IT), + AT) > [+2] =(1 + AT/T,) > 
AT/T,= (1 + ¿ui 
AT,=[(1 +8)" - 1]7, = (1/4)€'T,. (20) 


La Figura 8 muestra gráficamente la ec.(20). Note que un aumento, 
por cualquier causa, del flujo IR de un 0,01 del flujo inicial (o sea, de 
un 1%), produce un aumento en la temperatura de equilibrio radiativo 
de la Tierra de aproximadamente 0,7 K, o casi de 1 K. 
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FIGURA 8. Cambio de la temperatura efectiva terrestre debido a una variación de la 
alteración (fracción) del flujo ir saliente inicial. 


: Una de las cuestiones que es conveniente aclarar a estas 
alturas, es la interpretación que hay que darle a estos cambios de 
temperatura del orden de 1 K. Se podría argumentar que un cambio 
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de un grado Kelvin, es un cambio muy bajo como para que produzca 
cambios en el sistema climático. Vale decir que esta variación de 
temperatura es menor que aquella que ocurre normalmente entre el 
día y la noche sin que esto produzca trastornos climáticos. El punto 
es que cuando se habla de cambios en la temperatura de radiación 
(efectiva, atmosférica y superficial), se está hablando, respectiva- 
mente, de una temperatura promedio anual global (“temperatura 
climática”). Y una variación pequeña de estas temperaturas 
climáticas significaría variaciones mayores en las temperaturas 
(meteorológicas) que normalmente se presentan diariamente y que 
se usan para hacer los promedios globales. Apoyémonos en una 
analogía para entender esto mejor. 


Considere, por ejemplo, la nota definitiva, ND, obtenida por un 
estudiante, como resultado de promediar sus notas parciales. 
Supongamos que fueron cuatro exámenes parciales (EP, ¡= 1, 2,3 y 
4) que éste presentó. Así, ND = (EP, + EP, + EP, + EP,)/4. La nota 
definitiva crítica es ND = 09 ya que se puede esperar que se solicite 
un aumento de un punto para llevarla a diez y aprobar. ¿Qué 
significa elevar en un punto el valor de ND? Despejemos de la 
ecuación anterior cualesquiera de las notas parciales, por ejemplo, 
EP, = 4ND - EP, - EP, - EP,. Hagamos un ejemplo numérico: 
pongamos que para que el estudiante le dé ND = 09, las notas 
parciales fueron EP, = EP, = EP, = 10 y EP, = 06. Si se le aumenta 
en un punto a ND para que este valga 10, EP4 = 4x10 — 3x10 = 10 
puntos. En otras palabras, aumentarle de 09 a 10 la ND (diferencia 
de un punto), significa elevarle la nota del cuarto parcial de 06 puntos 
a 10 puntos, o sea, una diferencia de cuatro puntos, que puede 
equivaler a un problema que hizo mal o que no hizo. 


En relación a la aproximación (1), un modelo un poco más 
elaborado pero aun todavía sobre-simplificado, que considera la atmósfera 
dividida en varias capas, pero con e, = 1 para una más rápida captación del 


24 Introducción a la Física de la Contaminación Atmosférica y del Cambio Climático 


modelo (radiación de cuerpo negro), permite tener una idea un poco más 
realista del problema. 
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EL EFECTO DE INVERNADERO . 
EN UNA A ATMÓSFERA ESTRATIFICADA: 


La observación del perfil de temperatura troposférica, entre la 
superficie y la tropopausa, indica una disminución de ésta con la: 
altura en una cantidad que- oscila entre los 60 2C.y los 80 *C. Esta. 
disminución de temperatura puede ser explicada en principio con 
base a la transferencia de radiación en una atmósfera hipotética, 
discretamente estructurada en capas homogéneas plano-paralelas, : 
teniendo cada una asociada una temperatura T, (con ¡= 1,2... n), tal 
que T,> T, > T,> T, > ...> T,(Goody 8 Walker, 1975). 


Cada una de estas capas debe tener un espesor tal que la 
radiación emitida por una, debe ser absorbida por la adyacente. Así, 
las capas no deben ser ni demasiado gruesas ni demasiado 
delgadas. Dentro este contexto, una capa será demasiado gruesa si 
la radiación es emitida y re-absorbida dentro de la misma capa. Por 
el contrario, una capa será demasiado delgada si la radiación la 
atraviesa, sin experimentar ninguna absorción. Por lo tanto, cada 
capa tendrá un justo espesor necesario para absorber la radiación 
que le llega. Al aumentar la altura, la densidad de los gases 
absorbentes es menor, por lo que se espera que el espesor de las 
capas aumente también con la altura, permitiendo a la radiación 
viajar más lejos antes de ser absorbida. 


El número total de capas en. que una atmósfera puede ser 
dividida de esta manera, o bajo este criterio, se llama la opacidad o * 
espesor óptico, 1, de la atmósfera. Consecuencialmente, este 
espesor óptico depende de la cantidad de atmósfera que haya y de la- 
eficiencia con que los gases atmosféricos absorben la radiación IR. 
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En el equilibrio consideremos el balance de energía entre la 
diferentes capas presentadas en la Figura 9 (Goody €: Walker, 1975). 


| 1 | 
oT, Capa 2 (T,) 
! y 
a oT, Capa 3 (T,) 
E 
> CE — al EN z Capa 4 (7 ) 
Ao 
oy pa 
SUPERFICIE (T) 


FIGURA 9 
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La capa 1, radiativamente hablando hace satagio” con .el 
Diogo BxlScian Lal lo que: 


OT = aT,, en 
de donde, 7, =T, | E 
El equilibrio djativo: en la capa 1 sugiere también que: 
do TA + PT. = No Tos > T= 2T,= = 272, (22. 
lali si se establece el balance de energía en la capa 2: 
ON Ms toTA4 e => sa = -2T,- TA. (23) 
Haciendo uso de la relación (22) y (21), 
| Té ATA TS TARTA ea 
En la capa 3 el dins e de energía es: as 
Jo T, + oT,- Yo T¿+ doT, > T/= ato. r; (25 
Aplicando simultáneamente las relaciones (24) y 22). 
TE = 67% er; => 7, = ar. eS 
A comticudción en a capa 4: 
doT+ toT; toros lor, > > TA=2T2- TA en 
Por. aplicación de (26). y (24), 


T =8T;- ari T,= BT, (28) 
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Finalmente, en la capa 5: 
doT¿4ToT,=PoT¿+ do T¡¿= Ti =2T;- T,.. (29). 
Con la relación (28) y (26), la ec. (29) se transforma en, : 
7! =10T/- 475 T*= 67. (30) 


Las expresiones (22), (24), (26), (28) y (30) nos permiten encontrar 
una ecuación general de la temperatura superficial, 7, en función del 
espesor óptico en el infrarrojo, 7,,, de la atmóstera: le 


T.= (1 + número de capas)T, => T,=(1+1,)"T,, (31) 


La variable 1,,, como se deduce de la expresión anterior, es 
adimensional,,es decir, no tiene unidades. En el ejemplo que nos 
ocupa, el número de capas es igual a 5, por lo que el espesor óptico 
en el IR es igual a 5, con cada capa teniendo un espesor óptico 
individual igual a uno. Note, sin embargo, que a pesar de esto el 
espesor geométrico (altura) de la capa más profunda, la que está en 
contacto con. el suelo, es el menor puesto que esta capa es la más 
densa. A medida que subimos de capa, la densidad disminuye por lo 
que para que tenga igual opacidad óptica, el espesor geométrico 
debe aumentar. La disminución de densidad debe ser compensada 
por un aymento del espesor geométrico de la capa; de allí que en la 
Fig. 9 la.capa 1 es la más gruesa y la capa 5 la más delgada. Pero 
todas, individualmente, tienen una opacidad igual a uno. Por lo 
anterior, se puede afirmar que, en general, el espesor óptico es 
proporcional a: () el espesor geométrico x de la capa (o altura en 
este caso), (íf).la densidad p de los gases contenido en ella, y (ii a la 
“eficiencia” k con que estos gases extingan (absorción + dispersión) 
energía radiante; por consiguiente, r = kpx. Como esta variable no 


MEA 


tiene dimensiones, y considerando que la unidad de densidad puede *: 


ser expresada en [g/cm'], y el espesor plc en visi la unidad 
de k puede venir dada en [cm”/g]. 


Es importante entender, entonces, que por esta última - > 


definición, «la: opacidad o espesor óptico es un número que no: “* 


necesariamente tiene que ser entero (como el valor de:5 de'la Fig. 9). 
Considerando que T, = 255 K y que 7, = 288 K, el cálculo muestra 
según la. ec. (31), que la opacidad media de la atmósfera en el 
infrarrojo es de aproximadamente 0,63. Evidentemente un “aumento 
de esta opacidad por encima de este valor, debido a un aumento 
antropogénico de estos gases (por ejemplo, CO,), conllevará a un 
aumento de T,. ¿Habrá aumentado o estará aumentado sostenida- 
mente la opacidad de la atmósfera en el IR? 


La. ec.(31), al igual que la ec.(10),' son cada una 
representación matemática del efecto de invernadero natural de la 
atmósfera terrestre, en función respectiva de el espesor óptico en el 
IR (r,) y la emisividad (s,). Si la Tierra no tuviera atmósfera, o 'no 
contuviera gases absorbentes en el infrarrojo (gases invernadero), 1,, 
= e, = 0, y de'ambas ecuaciones, 7, = T,: no habría calentamiento 
global natural de la superficie. Pero como esto no es cierto, 1,, > 1, £, 
> 1 por lo que 7, > T,. En este último caso, una alteración de esta 
condición (recalentamiento A variación de la temperatura de 
radiación del globo terráqueo, 7,, como cuerpo en el espacio) puede 
venir expresada por la ec. (20) la cual depende precisamente de e (= 
ASS, que, como se ha indicado, representa la fracción dela 
variación de S, en función del flujo original), debido a una alteración 
en la concentración y/o composición del aire. Si no hay tal alteración, 
AS, = = 0, e =0 y por consiguiente, AT, = 0. 


Según la ec. (10), si la atmósfera fuera considerada como 
constituida por una sola capa (Fig. 5) y emitiera como cuerpo negro 
(e, = 1) > 7, = 1,1897, (o T/T, = 1,189). Por la ec.(31), para una sola 
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«capa (1, = 1), T/T, = [1 + 1]% = [2] = 1,189. De momento este 

cociente es T/T, = 288 K/255 K = 1,129 lo cual quiere decir que la 
atmósfera ni emite como cuerpo negro ni debe tener una opacidad 
óptica infrarroja de 1 (sino de -0,63 como se ha deducido arriba). 
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EL EFECTO DE INVERNADERO EN UNA ATMÓSFERA 
DE UNA SOLA CAPA CON MÚLTIPLE REFLEXIÓN 


tl dE 


Con atención. a las aproximaciones a, (6). y (8) lechas: anterior- 
mente, podemos considerar la posibilidad de que suceda lo contrario 
a lo indicado por ellas: (a) no toda la energía solar. llega a la 
superficie, parte de ella es absorbida antes y parte reflejada; (b)-la 
radiación no es el único. mecanismo de' transporte de energía; (c) no 
toda'la: energía que llega a'la superficie es absorbida por ésta. Cón: 
relación a estos literales definamos, respectivamente, por fla fracción 
de la energía incidente en. el visual no-reflejáda que es absorbida por 
la atmóstera (abosrtividad); por S,, el flujo no-radiativo (convección y 
calor latente de evaporación) proveniente de ' la: superficie; y por r, el 
albedo de la superficie [luego, (1 — r,) es la parte que absorbe la 
superficie y r, es la parte que refleja]. A su vez la capa atmosférica 
posee su albedo, que podemos denominar por r, [luego, (1- r,) es la 
parte que ingresa a la atmósfera y r, la parte reflejada al espacio] La 
Figura 10 ilustra lo que vamos a explicar. 


S, 1S, 
y ? 
y p « 
(1 z r)S, => 1 y ms, ira. reflex.atmós. a Sra. refiex.atmós, 
y ES, > (1-5)r.S, LES, > (11, IS, ES > 
S,= (1 - N(1 - 1IS, -> +48 ; > Py Py : 


1ra. reflex. sup. 2da.reflex.sup. 3ra. reflex.sup. 


(1 A r)S, SP (1 = IIS, Ea (1 ES r)rir7S, E (1 5 rreris, 






FIGURA 10. Múltiple reflexión en una atmósfera de una sola capa. 
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Inicialmente S, (recordemos que vale 342 Wm”) se distribuye 
en r,S, que se refleja al espacio y en (1- r,)S, que penetra en la 
atmósfera; de esta última, f(1— r,)S, es absorbida por ella y (1- f(1— 
r,JS, = S, llega a la superficie. De esta última, (1-- r,)S, es absorbida 
por la: superficie y r.S, es devuelta a la atmósfera (primera reflexión 

superficial). De ésta, (- r,)r.S, es absorbida por la atmóstera y 1,S, 
es devuelta hacia abajo (primera reflexión atmosférica). Al llegar ésta 
a la superficie, (1- r,)r,r,S, es absorbida por ella, y r,r.S, es devuelta 
a la atmósfera (segunda reflexión superficial). De ella, (1-r,)r,rÍS, es 
absorbida y r/r/S, es devuelta hacia abajo (segunda reflexión 
atmosférica). De. esta cantidad, (1- r,)r?r3S, es absorbida por la 
superficie y rórSS, es devuelta a la atmósfera (tercera reflexión 

superficial). e esta última cantidad, (1- r,)r,7r?S, es absorbida por la 
atmósfera y r//r?S, es devuelta por ésta hacia abajo. La m reflexión 
superficial (hacia arriba), siguiendo las expresiones en amarillo de la 

Fig.10, tiene la:forma general r,”*r."S,(%), y (1— r,)r,”"1,7'S (Y) es la 

expresión general para la m-ésima absorción de la superficie. La n 
reflexión atmosférica (hacia abajo), siguiendo las expresiones en azul 

turquesa en la Fig.10, tiene la forma general r, "r2S (Y, y (1- rr 

'rS, (>) es el término general para la n-ésima absorción de la 
atmósfera. Después de m y n reflexiones, con m > «w y N > 00, 

tenemos que, S 


Absorción Total en la idas paa a- INEA) S, = (1 — rr) 
(1 - rJ8, : = 9 


= k,(1 198, = k,(1=1)/(1= (1 - 1,JS,= BS,, (32) 
donde k,, = (1—r, 1)" y B= k,(1— 1) (1 -A(1- m6) 
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Absorción Atmosférica por Reflexiones = 3”,_, (1 —r,)r.(r,1)"S, = 
(1 — ELL — JAS, = K, (1 —5,5.S, 


3. rr > NA - 18, = k, MES y r0- S.. (33) 
Absorción Total de la Atmósfera = LK (1 - 0 ki : p+ Ñ - JS, 
| AS, (84) 
donde. A= Ik, (1 —1y rn > pr +(1 - ral | | 
Con las expresiones (32) y (34), podemos construir, haciendo 
un balance de energía, un sistema de dos ecuaciones con dos 
incógnitas donde ellas son, respectivamente, 7, y T,. - 
El balance de energía en la atmósfera es: 
AS, + €, 07, + hS, = 2e,07,,, (35) 
donde HS, = S,,. El balance de energía en la superficie es: 
BS, + £,07, = 07, +hS,. (36) * | 
De la ec. (35), 
T/ = (1/26,0)(0,0T; + AS, + hS,). (87) 
Sustituyendo ésta en ec. (36), 
| 7 =(S/K(2 -eJo]) (28+ A - h). (88) 
Sustituyendo en ec.(38), o | 


T? = (Se, o) M[1- e,Jh + e,8+ AJ). (89) 


34 Introducción a la Física de la Contaminación Atmosférica y del Cambio Climático. 


Sif=0,€,=1,r,=1,¿=0y h=0, la atmósfera no absorbe: 
nada (transmite todo en el visual), emite como. cuerpo negro, ella y la 
superficie reflejan todo en el visual y absorben todo en el IR, y no hay 
convección ni calor latente de evaporación. En este caso, B=1 y A= 
0; en otras palabras, tenemos el ejemplo de la Fig. 3. Aplicando 
ecs.(38)-(39) con estos valores, S, = 07; ='S, y T, = 279 K; 28, = 684 
Wm? = o T¿ = S, y T, = 331 K, lo que significaría un "perfecto efecto 
dé invernadero”. Lo anterior nos indica que las expresiones 
matemáticas (38) y (39) para las temperaturas están bien. Si 
cambiáramos r, al valor de 0,35, S, = 222 Wm” y S, = 444 Wm”, 
como en la Fig.4, dando, entonces, 7, = 250 K y T,¿ = 298 K.,... 


Problema 3: 


Con las ecs. (38)-(39), demuestre que bajo las condiciones 
dadas al principio del párrafo anterior, T, = [1/(2)*]7, [ec. (8)]. 


Es interesante mostrar cómo, simplificando al máximo el 
modelo de reflexión múltiple, o sea, poniendo solo r, = r,= 0, la 
misma atmósfera puede actuar como un escudo térmico en lugar de 
ser un amplificador del calor (Knox, 199). Estamos ahora hablando 
del caso de reflectividad cero en el visual para ambos componentes 
del sistema (atmósfera y superficie). Quedan en juego, entonces, f y 
€, Para ello, veamos cómo varía la absortividad-emisividad en IR (e,) 
como función de la absortividad en el visual (f), para diferentes 
valores (constantes) del parámetro libre T, (isotermas). De la ec.(38), 


e, = 2-[8,(2 —f/0T/]. (40) 
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FIGURA 11. Selecto conjunto de isotermas (Ts) calculadas con la ec.(40) cuando los 
parénetros s, y f varán entre O y 1 (Knox, 1999). a 


Para T, = 279 K tenemos la diagonal isotérmica que va 
subiendo de izquierda a derecha, desde el punto (e,,f) = (0,0) al punto 
(1,1). Para el primer punto tenemos el ejemplo de la Fig.1 (no hay 
atmósfera). Pero si f = e,, todo el tiempo como en esa diagonal, T, 
siempre será igual a 279 K (máximo efecto de enfriamiento o de anti- 
invernadero); y es como si la atmósfera no existiera porque lo que 
absorbe en el visual lo devuelve en el IR aun absorbiendo todo en el 
visual y devolviendo todo en el IR [punto (1,1)]. Cuando nos 
desplazamos de (0,0) a (1,0) por el eje de las ordenadas, la 
atmósfera no absorbe nada en el visual pero absorbe y devuelve en 
el IR. Con e, = 1, tenemos la atmósfera como cuerpo negro y T, = 
331 K que es el caso de la Fig.3 (máxima 7, o efecto de 
invernadero). Cuando nos desplazamos del punto (1,1) al (0,1), la 
atmósfera se hace más transparente al IR, manteniendo máxima 
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opacidad en el visual; como resultado, 7, cae a su mínimo valor de. 
234 K (mínimo efecto de enfriamiento o de anti-invernadero, e, = 0), 
porque la superficie recibe menos energía IR. En resumen, si la 
_ atmósfera se hace más transparente al IR, para un f fijo, T, tiende a. 
disminuir de 331 K a 234 K y viceversa (línea vertical a trazos en la 
Fig.11). Si bajo la misma condición f varía, T, también tiene esa 
variación y viceversa (línea diagonal a trazos en la -Fig.11). Esto 
muestra claramente cómo la radiación IR juega un papel Amporante 
en el efecto de invernadero o calentamiento global. 


Problema 4 E | 
Ñ Grafique- e, VS f para diferentes isotermas atmosféricas, 


compárelas con las del gráfico de la Fig.11 y discuta los 
resultados. 


A 5: 


Grafique e, vs f para diferentes isotermas superficiales, con r, = 
0,25 y r, = 0,16 y discuta los resultados. 
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CONSECUENCIAS. DEL RECALENTAMIENTO GLOBAL: 
PARÁ MUESTRA, UN BOTÓN | 


El aumento sostenido de los. peeestara de la metan. dé Mob 
invernadero, como . el CO, (dióxido o: bióxido .de. carbono), CH, 
(metano), NO, (dióxido de nitrógeno), etc., por acción antropogénica, 
conduce, como se ha ilustrado, a un aumento de las temperaturas 
efectiva del planeta, atmosférica y superficial, respectivamente, como 
respuesta..climática del sistema atmósfera-superficie mistno. .Estas 
respuestas climáticas tarda su tiempo en manifestarse. En: el caso 
del recalentamiento global, este puede tardar décadas: en aparecer 
(Tabla 1). Este aumento de concentración de gases invernadero :'ha 
venido ciertamente sucediendo desde que comenzó la era: industrial 
en el siglo XIX, por lo que los cambios climáticos que se ha venido 
observando, se le han atribuido a este.efecto. Y esto está siendo 
sometido a-estudio e intensos debates científicos. Para ello se hacen 
proyecciones a futuro, bajo diferentes escenarios y Baras, a 
base de modelos matemáticos. :- “4 da 4 


En todo caso, este aumento de temperáturas trae como 
consecuencia la potencial alteración de' prácticamente todo'lo que 
ocurre en la atmósfera, litósfera, hidrósfera; biósfera y «criósfera; 
ejemplo de ello sería ta modificación del ciclo hidrológico. (lluvias, 
nubosidad, sequía, etc.), aumento del nivel del. mar, descongelación 
de los polos y otras cosas más (Erickson, 1992). 


- Tomemos como ejemplo. la descongelación de los casqueles 
polares y el aumento del nivel del: mar. Mostraremos con un sencillo 
ejemplo" cuál podría el orden. de magnitud «de este efecto como 
secuela a largo.plazo del recalentamiento global (Hafemeister, 1982). 


Si se supone que un aumento global de temperatura media, 
causado por un aumento artificial de la concentración de CO,, es 
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millones de años. Por esta razón, la variación del nivel de los 
océanos, a causa de este fenómeno, que inundaría. mundialmente las 
tierras bajas, es una posibilidad diferente. >=“ E 
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PREVISIONES Y PRECAUCIONES 


Hay ciertas evidencias que parecieran indicar y advertir que el 
recalentamiento global terrestre ya está en marcha. Después de 
haberse detectado grandes fracturas y grietas en la capa de hielo del 
Polo Sur y la creación, por ello, del iceberg más grande que jamás se 
haya observado, se ha llegado a sospechar que esto está siendo 
“provocado precisamente por un aumento del efecto de invernadero 
(este singular hecho de! iceberg, fue reseñado por los medios de 
comunicación social internacionales). Por otra parte, en la oscilación 
normal que tienen las temperaturas a lo largo de un año, se ha 
podido notar que en general las mínimas ya son tan mínimas, sino 
que han auméntado; en otras palabras, las temperaturas mínimas ya 
no están bajando a sus valores acostumbrados. Y esto también se le 
está atribuyendo al recalentamiento global. En todo caso, lo que sí es 
ciorio es que la concentración de gases invernadero en el aire por 
acción de la actividad humana (quema de combustible fósil, quema 
de biomasa, etc.) ha ido en aumento, y esto plantea, como se ha 
ilustrado en este trabajo, un potencial efecto sobre el clima. 


De ahí que una forma lógica de prevenir y atenuar o mitigar el 
problema es lograr una reducción de las emisiones de estos gases a 
la atmósfera. Para ello a escala internacional, a través de la 
Organización de las Naciones Unidas (ONU), los países miembros 
firmaron el 13 de ¿unio de 1992 un instrumento jurídico ambiental 
denominado “Convenio Marco de las Naciones Unidas sobre el 
Cambio Clirziático” cuyo objetivo general es (MRE/MARN, 2000), 


“Lograr, de conformidad con las disposiciones pertinentes de 
la Convención, la estabilización de las concentraciones de gases de 
efecto invernadero en la atmósfera a un nivel que impida 


MARCOS PEÑALOZA 41 


interferencia antropógenas peligrosas en el sistema climático. Ese 
nivel debería lograrse en un plazo suficiente para permitir que' los 
ecosistemas se adapien naturalmente al cambio climático. Asegurar 
que la producción de alimentos no se vea amenazada y permitir que 
el desarrollo económico prosiga de manera sostenible”; o 


Este convenio, que tiene grandes implicaciones económicas y 
sociales para nuestro país, .en virtud de su condición de país cuya 
economía es altamente dependiente de la exportación de combus- 
tible fósil, fue ratificado por Venezuela el 27 de Aer de 1994 
(Gaceta Oficial N* 4825 Extraordinaria). ' 


Para implementar este convenio, el mismo admite en su 
Artículo 17, la adopción de protocolos, y en su Artículo 15, la de 
enmiendas. En este sentido, actualmente está bajo intensa discusión, 
polémica el primero de ellos conocido con el nombre de “Protocolo 
de Kyoto” (MRE/MARN, 2000). Su objetivo general es: 


“Comprometer a los países industrializados a la reducción de 
sus emisiones de gases de efecto invernadero por lo menos en un 
5% con relación a los niveles de 1990, para el período entre el 2008- 
2012. 


- Estos dos instrumentos, junto con el “Convenio de Viena para 
la Protección de la Capa de Ozono” (MRE/MARN, 2000), firmado el 
22 de marzo de 1985 (ratificado por Venezuela el 19 de julio de 1988; 
gaceta Oficial N2 34010), y su “Protocolo de Montreal Relativo a las 
Sustancias Agotadoras de la Capa de Ozono” (MRE/MARN, 2000), 
firmado el 16 de septiembre de 1987 (ratificado por Venezuela el 11 
de enero de 1989; Gaceta Oficial N* 34134), constituyen hasta el 
momento los mayores esfuerzos hechos a escala mundial para 
detener el daño que se le está haciendo al ambiente atmosférico, 
para protegerlo y para conservarlo. 
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El conocimiento detallado de estos documentos permite tener 

una idea completa de estos esfuerzos y de los obstáculos políticos, 
económicos, sociales, etc., que se han superado y que todavía hay 
que superar, para aplicarlos. Esto se deja como actividad opcional y 
complementaria al lector; hecho esto, sería interesante deducir, 
- entonces, cuál es la diferencia entre un convenio y un protocolo... . 
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CONTAMINACIÓN DE UNA CIUDAD 
(MODELO DE CAJA FIJA) 


- La calidad del aire de un cierto ambiente es un aspecto de la 
contaminación atmosférica, que debe ser tratado científicamente 
para Colaborar en la toma de decisiones políticas y de planificación 
respecto a la gestión y desarrollo de infraestructuras e industrias. 
Como análisis de, sistemas, la calidad del aire puede ser 
representada de acuerdo a la Figura 12 (Kiely, 1999). 


, Calidad del Aire 
———>P| Modelo Matemático Previstas 


- Meteorología |- 


- FIGURA 12. Análisis del sistema para estimar la calidad del aire. 


De los modelos matemáticos, el modelo de caja fija es uno de 
los más simples que se utilizan para hallar una primera estimación de 
los valores de la concentración de los contaminantes. Se basa en la 
conservación de la masa de un contaminante en un espacio físico del 
tipo caja. Considere una ciudad rectangular, como se muestra en la 
Figura 13, con el fin de calcular la concentración del contaminante 
del aire en esta ciudad (De Nerves, 1998). 
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Velocidad del viento 1 


FIGURA 13. Modelo de ciudad rectangular para aplicar el modelo de Caja Fija en la 
p - estimación de su contaminación. 


Para simplificar el análisis, hagamos la siguientes aproxima- 
ciones: 


(a) La ciudad es un rectángulo con dimensiones Wx L, y con uno de 
«Sus lados paralelos a la dirección del viento. 


(oj La turbulencia atmosférica produce el mezclado completo y total. . 
de. los contaminantes hasta la altura de mezclado H, y no hay 
A por encima de esta altura. 


(0) Esta. turbulencia es bastante fuerte en la dirección contra el 
viento, de modo que la concentración del contaminante es uniforme 
en el volumen completo de aire que está sobre la ciudad y no a más 
altura, en el lado en la dirección del viento y en el lado contra el 
viento. Esta hipótesis es bastante contraria a lo que se observa en la 
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naturaleza, pero permite una gran simplificación del . aspecto. 
matemático. : 


(d) El viento sopla en la dirección x con velocidad u. Esta velocidad : 
es constante'e inde- pendiente del tiempo, lugar o elevación por: 
encima del suelo. También esto es contrario a la observación; la 
velocidad del viento aumenta con la elevación. En lo que. sigue se 
usa la u promedio entre aquélla al nivel del suelo y la 
correspondiente a la altura H. : A od 


(e) La concentración del contaminante en el aire que entra a la 
ciudad (en x = 0) es constante e igual a b [b por concentración "de 
fondo" (background en inglés)]. En este modelo, como en la mayoría, 
las concentraciones se dan en [g m”] o en [ug m?]. 


(f) La tasa (índice) de emisiones del contaminante del aire de la 
ciudad es Q [g s*]. Normalmente éste se da como un índice de 
emisiones por unidad de área q [g s' m”]. La expresión Q = gA 
permite convertir uno en el otro, donde A es el área de la ciudad (= 
WL). Este índice es constante y no cambia en el tiempo. 


(g) Ningún contaminante sale o entra por la parte superior de la caja, 
ni por los lados que no sean paralelos a la dirección del viento. 


(h) El contaminante tiene una duración suficiente en e aire como 
para Au5 su velocidad de destrucción sea cero. 


En este modelo, como en casi la mayoría, el balance de 
materiales es la idea que subyace en la física de los mismos (ley de 
“conservación de la masa). Se entiende como un. balance "de 
materiales, una especie de contabilidad en la que se aplica una 
ecuación general del balance a alguna especie contaminante. Dadas 
ciertas condiciones de frontera o contorno, esta ecuación se pad 
escribir como (De Nerves, 1998), 
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(Velocidad de Acumulación) = (Todos los Gastos de Entrada) — 
[Todos los Gastos de Salida) + ielecidad de Creación) + (Velocidad 
de Destrucción). (41) 


Esta ecuación, bajo las restricciones arriba impuestas, se reduce a: 


0= [Todos los Gastos de Entrada) - [Todos los Gastos de Salida) > 
(Todos los Gastos de Entrada) = (Todos los Gastos de Salida). (42) 


Un gasto de entrada del contaminante hacia dentro de la cara 
contra el viento es 


Gasto de Entrada 1 = uWHb. (43) 


El factor uWH expresa el volumen de aire que cruza la frontera 
“contra el viento del sistema por unidad de tiempo; note que esta 
- expresión tiene unidades de [volumen/tiempo] que al BEAR 
por la concentración de fondo da [masa/tiempo]. 


El otro gasto de entrada es aquel proveniente de la propia 
ciudad, que entra al aire por la frontera, o cara, inferior de la caja; 
este es igual a . 

Gasto de Entrada 2 = Q = qWL. (44) 
La suma de la ec. (43) y la ec. (44) da todos los gastos de entrada: 


(Todos los Gastos de Entrada) = (uHb'+ qL)W. (45) 


Con relación a los gastos de salida el único que sale es aquel 
- que lo hace hacia fuera por la cara del extremo derecho a la caja (en 
x =L), perpendicular al viento. Como en el gasto de entrada 1, el 
gasto de salida es el volumen uWH pero multiplicado... por la 
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concentración c (en L): . 
[Gasto de Salida) = uWHc. (46). 


Sustituyendo ec. (46). y ec.(45) en, ec. (42),. y. despejando..c,..se 
encuentra que 


c=b+ qluH. (47). 


Es fácil ver que la ec. (47), como se muestra en la Figura 14, 
representa una recta donde la variable dependiente es c, la variable 
independiente es q, el término: (L/uH) -es la peniligote;! yb: es el 
término independiente. , 


Ñ 
| 
| 
| 





Concentración - 
en cl aire ambiemo 


Concentiación de fondo, + 


AAA 


A 
Índice de emisiones, q, g/s kay? 0 


FIGURA 14. Representación gráfica de la ecuáción (47) (De Nerves, 1998). 


El modelo predice que la concentración en el aire es 
directamente proporcional al Índice de emisiones q de la, ciudad; a 
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mayor emisión, mayor contaminación. Pero también nos dice con que * 
rapidez esta contaminación. se presenta; y esta viene dada por la 
pendiente la cual, a Su vez, es directamente proporcional al largo L 
de la ciudad e inversamente proporcional a la velocidad del vie nto 17 y 
a la altura H de la capa de mezclado. Si L aumenta, la rapidez de la 
contaminación sube. Si u aumenta, la rapidez baja ya que el viento 
limpia más eficientemente el aire. Si H aumenta, los contaminantes 
están más diluidos. por, lo que | la rapidez baja. 





Pre blema 6. 


(a) Dejando q y lo demás constante, muestre escquemáticamente 
con un gráfico, cómo varía c con el largo L' de la ciudad y 
discuta cualitativamente el resultado. (b) Dejando q y lo demás 
constante, -muestre esquemáticamente con un gráfico, cómo 
varía c con la velocidad' u del viento y discuta cualitativamente el 
resultado. * “Dejando q y lo demás constante, muestre 
esquemáticamente con un gráfico, cómo varía c con la altura H 
de la capa de mezclado y discuta cualitativamente el resultado. 


El modelo-en cuestión es una gran simplificación de lo que en 
la realidad:sucede, pero muestra con sencillez cómo influyen en una 
primera aproximación las variables más importántes. Note cómo el 
modelo no incluye a W como variable y esto se explica porque si se 
duplica: el “ancho W el la ciudad; ilegal de coristante, * cono 
cambia. * A nera 


Que las concentraciones en los bordes de la ciudad, en la 







peor. de. las. simplificaciones. Le sigue la á proximación (b) puesto que 
no todas las. emisiones son las mismas a lo largo. y ancho de la 
ciudad; pero si el mezclado de los contaminantes en el aire nó es 
perfecto, ésta se vuelve la peor. Por ejemplo, para mayor parte de los 
contaminantes, q puede ser bajo en los suburbios y mucho más alto 
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en las zonas industriales, pero el buen mezclado hace que.asumir q 
, constante, en toda la ciudad sea una buena suposición. 


Hagamos un ejemplo numérico aplicando el modélo de caja 
fija: Una ciudad y su atmósfera tiene las siguientes características: W 
=5 km, L=15km,u=3ms' y H= 1000 m. La concentración contra 
el viento, o de fondo, de, por ejemplo, CO es b=5ugm” y la emisión 
de contaminantes por unidad de área de la ciudad es q = 4x10* ( gs” 
m?. ¿Cuál es la concentración dé (CO sobre la ciudad? Apicándo la 
ec. (47) : 


c=5 gm” + [4x10* g s* m”][15000 mA3 m.s” x 1000 m)] =(5 +20) 19 m? =25 Hg m. 


Vale la pena hacer resaltár que manteniendo. lá velocidad del 
viento constante, nos estamos refiriendo” a una. condición” -meteo- 
rológica específica. Si quisiéramos referimos a diferentes concentra- 
cioniés para diversas circunstancias, a lo largo de un año, por ejem- 
plo, tendríamos que estimar un promedió, ánual de la concentración 
de algún contaminante con base a la distribución de frecuencias de 
diversos valores de u y H, calcular cada concentración respectiva con 
la ec.(47) y, a continuación, multiplicar por el número de veces que 
esta se presenta (frecuencia) y sumar, 0 sea, 


¿Concentración Promedio Anual =S Concentración pe Condición Metaorlógca) 
: Sobre todas - . 
las condiciones. 
meteorológicas a 
x (Frecuencia de Ocurrencia de esa Condición Meteoroiógica). (48) 


Hagamos un ejemplo numérico. Para el ejemplo anterior, las 
. condiciones meteorológicas descritas (u=3ms',H= 1000, m) 
ocurren el 40% de las veces. En el 60% restante, el viento sopla 
formando ángulos rectos con la dirección mostrada en la Fig.13,.con 
una velocidad de 6 m s* y la misma altura de mezclado. ¿Cuál es la 
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concentración promedio anual de CO en esta ciudad? Ya tenemos la 
concentración para el 40% de las veces (25 g m”). Calculemos la 
concentración para el 60% de las veces: 


Co = 5 ug m* + [4x10* g s* m*][5000 m/(6 m s'x1000 m)] = 8,83 ¡ig m” . 


Obsérvese que la dirección del viento ha intercambiado los valores 
de W y L. Si se usa este último valor, más el del ejemplo anterior, 


(Concentración Promedio Anual) = (25 jg m*x0,4) + (8,33 1g m*x0,6) = 15 yg m* 


Una: aplicación más realista de la ec.(48) requiere sumar 
sobre una gran cantidad de condiciones meteorológicas y sus 
correspondientes Índices de emisiones, en- lugar de: las.dos del 
ejemplo anterior. En ciudades con un gran consumo de calefacción 
doméstica, las fuentes de partículas producen Índices mayores de 
emisiones poriunidad de área en el invierno que en el verano. Los 
índices de emisiones de otros contaminantes pueden variar en 
lapsos de horas o de días. Las emisiones del parque automotor son 
muchos más altas en las horas picos que a cualquier otra hora del 
día. Por estas razones, la ec.(48) tendrá que modificarse para tomar 
en cuenta estas variaciones en las emisiones. 


Una. de las situaciones de mayor interés que se desearía 
conocer, es aquella relacionada con la concentración más elevada. 
Es necesario tener información, entonces, de la velocidad del viento, 
su dirección y sentido, la altura de mezclado y la concentración de 
fondo, que correspondan al peor caso. Y no siempre esta información 
está disponible. Una manera de hacer esto, bajo esta condición, es la 
propuesta principalmente por Larsen (1961) a: través de una 
manipulación sencilla de la ec.(47). Manteniendo constante a u, Ly H 
regresemos a la Fig. 14. Si deseamos reducir la. concentración de c, 
a C, es posible calcular con facilidad el índice permisible de 
emisiones q como 
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“a q (4 UH! Lo, 49) 


Si < se conocen 9. yo , Se despela. la. pendionte UH. (a, - pa 
ialacós en ec. 2.49) se ou0ane e, : 


q/q,= (6,- — bic, —b). 60). ña 





Note que los parámetros meteorológicos ya no pee porque se 
han eliminado, bajo el supuesto que estos son los. mismos para la 
-situación presente (c,,q,) y: la situación futura (c le 2%). Apoyados en, la 
ec. (50), la reducción fraccionaria. que $ se. necesitaría en el índice 
viene dada por : 


(Reducción. Fraccionaria en el Indice de Emisiones) = (q aa = =, 1 - ala 








1-10, MG 09606 SN, 





Hagamos un ejemplo numérico. En un país desarróliódo; la 
norma de calidad del aire ambiente para las partículas en 
suspensión, en 1971, era en promedio anual dé 75 ¡yg m”: “En 1970, 
la concentración promedio anual de partículas, medida, en una 
estación de monitoreo en el centro de una gran ciudad de ese país, 
fue de 190 ug m”. La concentración de fondo” se estimó en 20 Hg m?. 
¿En qué porcentaje tendría que reducirse el Índice de emisiones de 
partículas, por debajo del nivel de 1970, para cumplir con la norma 
de calidad del aire ambiente de 1971? 


IET Fraccionaria Sar el Indice de Emisiones) = Al 90 - 7500 - 20)= -0, 67, 


es decir, en un 67%. TA ec. (51) se conoce como el iodo 
proporcional o la ecuación de repetición y se ha usado ampliamente 
para calcular la reducción necesaria en los índices de emisiones, en 
procura de cumplir con las normas de calidad ambiental en lo que al 
- aire se refiere. Tiene la ventaja que es sencilla y normalmente 
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requiere datos de entrada disponibles con facilidad. Por el contrario, 
debido a su gran simplificación no se espera que dé predicciones 
exactas. De acuerdo al ejemplo anterior, los porcentajes calculados 
de reducción con esta ecuación tienden a grandes, lo que no hacen 
otros modelos más complejos. En consecuencia, este modelo se 
aplica con mayor frecuencia a ciudades con problemas de 
contaminación menos graves que aquellas con mayores. 


Un inconveniente de los modelos de caja fija es que no hacen 
«distinción entre grandes números de pequeñas fuentes que emiten 
- sus contaminantes a baja altura (como los carros, las casas, los 
basureros, etc.), llamadas fuentes de áreas, y el pequeño número de 
grandes fuentes que emiten cantidades mayores a mayores alturas 
(chimeneas industriales, etc.). En una simplificación sencilla, para 
hallar q ambos tipos de fuentes se suman. Pero en un enfoque más 
real, hay que tomar en cuenta que a mayor altura de la fuente, la 
concentración al nivel del suelo disminuye aunque más lejos de ésta 
la concentración aumente, y estos detalles no pueden ser tratados 
únicamente con un simple modelo de caja. 


Para una ciudad existente, que ha llegado a sus límites de 
crecimiento, nada se puede hacer para cambiar L, H y u. Pero si se 
está planeando una ciudad a futuro, el modelo de caja fija anterior 
_puede ayudar a comprender cómo se pueden manejar las variables 
- que intervienen en él, para minimizar la contafninación atmosférica 

que de ella se origine. Por ejemplo, el lado más ancho debe ser 
perpendicular a la dirección del viento (Wtan grande como se pueda) 
y el lado en la dirección al viento tan corto como se pueda (L tan 
.pequeña como se pueda). Si esto no se puede hacer, debiera 
escogerse un lugar donde H y u sean grandes. Lo anterior sugiere 
que las ciudades ubicadas en valles, perpendiculares a las corrientes 
de aire son las menos favorecidas por este criterio (pero las hay 
- bastantes; ¿puede dar ejemplos?). Sin embargo, las variables b y q 





:54 Introducción a la Física de la Contaminación Atmosférica y del Cambio Climático 


Sustituyendo, 
D= (4 14.000" x 500 ln / 1416 0.010 kl x 26h 900 si) 
¡Dz 3,590 km = 3.590 m1 


Este resultado es más de seis veces la altura de la estructura 
más alta construida por el hombre hasta la fecha y está más allá de 
- las capacidades actuales de la ingeniería estructural. Además, los 
“cálculos demuestran que la potencia requerida para impulsar el 
caudal sobrepasa la cantidad de energía eléctrica que se genera 
para. mantener una ciudad de este. tamaño. Como conclusión, el 
“problema de la contaminación del aire debe ser tratado. altemativa- 


mente bajo el procedimiento de la reducción de las emisiones. No 
* hay otra salida. 


Para cerrar esta sección, nos vamos a referir brevemente al 
. caso en que en la ec.(42) el balance de materia no sea cero. Esto 
significa que hay un aumento o disminución de la contratación con el 
tiempo (se rompe el equilibrio o balance), dependiendo del término 
. de entrada o de salida que predomine; en este caso, 


(da/dh) = [q/H] + (u/L(b — c). (52) 


La solución de esta ecuación corresponde a la de 'una 
ecuación diferencial de primer orden y de primer grado que aquí no 


corresponde resolver. Sin embargo, a continuación damos dicha 
solución: 


di 


ct) = c,_¿exp(utiL) + [(qL/uH) + BI — exp(-utL)]. (53) 


Esta ecuación muestra que si f > o, exp(-utlL) > 0 y c > [(qL/uH) + b] 


(modelo de caja fija); en t = 0, exp(utiL) = 1 y C = Cho 


MARCOS PEÑALOZA 55 


COMENTARIOS FINALES 


Con el entendido de que el problema es mucho más complejo, 
hemos hecho un intento para ilustrar y explicar tan simple como nos 
ha sido posible (ver aproximaciones hechas arriba), el fenómeno del 
recalentamiento global terrestre a través de un posible aumento del 
efecto de invernadero por contaminación artificial del aire, concebido 
básicamente como un aumento de la temperatura superficial del 
planeta con relación a su temperatura efectiva de radiación hacia el 
espacio. Lo hemos hecho apoyándonos en un simple modelo, que 
toma en cuenta la ley de la conservación de la energía y la ley de 
Stefan-Boltzmamn, y que considera a la atmósfera formada por una 
sola capa y por varias, en las que el balance de radiación entre la 
energía solar entrante y la energía infrarroja emitida por el sistema 
atmósfera-superficie, se ha establecido. Brevemente hemos revisado 
una de sus consecuencias por medio de un ejemplo sencillo, y 
hemos citado las gestiones que se están llevando a cabo para mitigar 
sus efectos. 


En segundo lugar hemos mostrado, a pesar de que también 
el problema es más complejo, cómo la contaminación atmosférica de 
una ciudad puede ser tratada desde punto de vista científico para 
establecer criterios de calidad de aire y ayudar a la toma de 
decisiones y diseños de política ambiental relacionadas con el medio 
ambiente. Esto se ha hecho en base a un modelo simple de balance 
de materiales, que toma en cuenta la ley de la conservación de la 
masa. 


Con estos aportes, creemos estar en consonancia y estar : 
honrando el compromiso adquirido por Venezuela, al ratificar el 
Convenio Marco de las Naciones Unidas sobre el Cambio Climático 
(MRE/MARN, 2000), que se encuentra contemplado en el Artículo 6 
(Educación, Formación y Sensibilización de! Público): 
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“Al llevar a la práctica los compromisos a que se refiere el inciso i) del 
párrafo 1 del Artículo 4, las Partes: 

a) Promoverán y facilitarán, en el plano nacional y, según proceda, 
en los planos subregional y regional, de conformidad con las leyes y 
reglamentos nacionales y según su capacidad respectiva: 

í. La elaboración y aplicación de programas de educación y 
sensibilización del público sobre el cambio climático; 

li. El acceso del público a la información sobre el cambio climático y 
sus efectos; 

li. La participación del público en el estudio del cambio climático y 
sus efectos y en la elaboración de las respuestas adecuadas; y 

iv. La formación de personal científico, técnico y directivo. 

b) Cooperarán, en el plano internacional, y, según proceda, por 
intermedio de organismos existentes, en la actividades siguientes, y 
las promoverán: 

i) La:preparación y el intercambio de material educativo y material 
destinado a sensibilizar al público sobre el cambio climático y sus 
efectos; y 

li) la elaboración y aplicación de programas de educación y 
formación, incluido el fortalecimiento de las instituciones nacionales y 
- el intercambio o la adscripción de personal encargado de formar 
expertos en esta esfera, en particular para países en desarrollo”. 
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